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We develop a consistent treatment of disorder effects in strong coupling superconductors. We 
use two different approaches, starting either from above or below the transition temperature, and 
show their equivalence. The normal state approach is a natural extension of a recently developed 
non-perturbative resummation technique. For the superconducting state we extend the standard 
Eliashberg theory to include disorder effects. We obtain a compact set of equations which are 
physically transparent and easy to use. They correctly reproduce the usual strong coupling results 
and disordered perturbation theory results in the appropriate limits. We numerically solve these 
equations to obtain transition temperature as a function of disorder for different phonon couplings. 
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The phenomenon of the suppression of the transi- 
tion temperature in thin superconducting films by disor- 
der has been extensively investigated in recent years [jj] . 
The physical mechanism is easy to understand: disorder 
causes electrons to move diffusively, leading to a less effi- 
cient screening of the Coulomb repulsion between them. 
This manifests itself as a low-momentum singularity in 
the dynamically screened Coulomb interaction. Naively 
one would expect this to lead to a very strong suppression 
of superconductivity due to the enhancement of the effec- 
tive Coulomb repulsion W. However this turns out not 
to be the case because of a cancellation between the vari- 
ous singular terms which arise, for example, in perturba- 
tion theory ||. This apparently fortuitous cancellation 
is in fact due to gauge invariance fjj, and the screened 
Coulomb potential can therefore be treated as if it were 
featureless. We note that all these subtleties are present 
even in the simplest possible calculation, namely first or- 
der perturbation theory assuming a featureless BCS at- 
traction. However this model is too restrictive for the 
following reasons: (i) perturbation theory is not valid in 
the limit of strong disorder often seen experimentally; (ii) 
many of the experimentally relevant materials such as Pb 
and Pb-Bi are strong coupling superconductors for which 
the details of the phonon spectrum are important [p|-[lO[ . 
It is therefore desirable to extend the simple theory to 
accomodate these features in a consistent manner. 

A method of going beyond perturbation theory has 
recently been developed by Oreg and Finkel'stein [pd] . 
These authors use a ladder summation technique which 
treats the BCS attraction and impurity-dressed Coulomb 
repulsion on an equal footing. This is possible because 
both interactions can be treated as featureless, as ex- 
plained above. The fact that this technique reproduces 
the renormalization group equation known from the in- 
teracting non-linear cr-model [ p"2| proves its essential cor- 
rectness. Its great advantage is its transparency and the 
fact that it is easy to use. The aim of this paper is 
to extend this technique to include strong coupling ef- 
fects. Initially one might think that this would further 
complicate the problem. Countcrintuitively the problem 



actually becomes simpler because one realises that the 
original technique is best understood in terms of strong 
coupling theory. Moreover this realisation allows us to 
naturally extend the non-perturbative treatment of dis- 
order to below the transition temperature. 

In this paper we present a consistent treatment of 
the effects of disorder and strong coupling on the tran- 
sition temperature in superconductors. We use two dif- 
ferent approaches which we show to be equivalent: (a) 
Starting from the normal state we derive the pair prop- 
agator and identify the transition temperature as the 
temperature at which it becomes singular; (b) follow- 
ing the Eliashberg theory of strong-coupling supercon- 
ductors, we start from below the transition temperature 
and we derive the superconducting self-energy and iden- 
tify the transition temperature as the point at which the 
order parameter goes to zero. The former approach has 
the advantage that it does not make any assumptions 
about the structure of the superconducting state and so 
is the natural way to initially address the problem. The 
latter approach has the advantage that it allows us to 
gain information about the superconducting state at all 
temperatures below the transition temperature at no ex- 
tra cost. The equivalence of the results obtained by the 
two approaches justifies the assumptions we make about 
the superconducting state. 

We begin by calculating the zero-momentum pair 
propagator, T(u>,u>'), shown in Fig. 1. This yields the 
equation 

r(w, uT) = r («, u')+tJ2 r (w, u")C(w")t(u", w') (i) 



where the uj are Fermi Matsubara frequencies, Tq is the 
elementary interaction vertex and C is the Cooperon im- 
purity ladder. All momentum dependences have been re- 
moved by averaging over the Fermi surface. We identify 
the transition temperature as the temperature at which 
r diverges. This will occur when there is an eigenvector 
4>(lo) of the homogeneous part of Eq. (|l|), which leads to 
the equation 
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If we define A(w) = T^^' u>')</)(w') the above equa- 
tion has the familiar Eliashberg form for the order pa- 
rameter. 



A(w) =T^r (w,w")CV')A(w")- 



(3) 



The elementary interaction vertex To contains the 
electron-phonon contribution, Tq^, and the Coulomb 
contribution. The latter naturally splits into two parts: 
the bare Coulomb interaction, To t c, for high momen- 
tum and frequency transfer, and the impurity-dressed 
Coulomb interaction, Tq^zs otherwise. The second con- 
tribution is important for low momentum transfers, \q\ < 
l/£, and low frequency transfers, \uj\ < 1/r, where I is the 
elastic mean free path and r is the elastic scattering time. 
This is because the disorder-screened Coulomb interac- 
tion keeps its long-range momentum singularity and its 
effective magnitude is therefore greatly enhanced. Gauge 
invariance leads to a cancellation of singular terms 
which allows us to treat the disorder-screened interaction 
as featureless with magnitude Vdis = 1 /2iV(0), where 
N(0) is the electronic density of states at the Fermi sur- 
face. 

The electron-phonon contribution, T$ph, is given by 
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where we have averaged over the Fermi surface in the 
standard manner. In the above equation <7k,k',A is the 
electron-phonon coupling for phonons of polarisation A, 
£>A(k — k',w — a/) is the phonon Green function and 
a 2 (z)F(z) = JV(0)£ A (| 3k ,^A| 2 <S(z r c. k - k <,A)) FS . The 
bare Coulomb contribution, Yq^Ci is just the Coulomb 
potential averaged over the Fermi surface, 

r 0tC = Vc = (V(k-k')) FS . (5) 

The disordered Coulomb contribution is given by 
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Fo,dis = Trilog 



(M + M)i 



(6) 



where we define the dimensionless measure of disorder, 
t = V dls /4ir 2 D = R a /R Q , and R = 2wh/e 2 = 162fcO is a 
quantum unit of resistance. 

First we include only the bare Cooperon ladder, 
Cq = nN(0)/\uj\, to yield the self-consistency equation 
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FIG. 1. Diagrammatic equation for the pair propagator 
F in terms of the elementary interaction vertex To and the 
screened Cooperon C. To consists of electron-phonon inter- 
action (zig-zag line), the bare Coulomb interaction (wiggly 
line) and the disorder-dressed Coulomb interaction block tA. 
£ is the self-energy correction to the Cooperon due to the 
impurity-screened Coulomb interaction. 
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where /i = N(G)Vc- The sums over the Matsubara fre- 
quencies in the above equation are cut off at the Fermi 
energy. However the phonon interaction is effective only 
for energies of the order the Debye frequency. We there- 
fore follow the standard procedure of introducing a 
frequency tu c above which only the bare Coulomb term 
contributes. The equation for A(u>) with u) > lo c then 
decouples and the only effect of these frequencies is to 
renormalize /i to 



[i = 



1 + fibx (ej?/o; c ) 



(8) 



We must therefore replace /j, by \i* in Eq. (^), and cut 
off the sums at frequency lu c . 

Next we consider self-energy corrections to the 
impurity-dressed electron Green function. Having av- 
eraged out the momentum dependence, the self-energy 
has the form = iu>[l — Z(u>)], The dressing of 

the electron Green function leads to the Cooperon being 
modified to C(u>) = ttN(0)/ (Z(ui)\ui\), as shown in Fig. 
(1). At this point we note the difference between the ef- 
fect of including the impurity self-energy in the original 
electron Green function, and the effects of including the 
phonon and Coulomb terms. The impurity correction, 
Si m p = {i/^ T )sgn(ui), is exactly cancelled in the process 
of constructing the impurity ladder - this is essentially 
the content of Anderson's theorem p3| which says that 
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non-magnetic impurities have no effect on superconduc- 
tivity at the mean field level. The contributions to the 
self-energy are 



(a) 



(b) 



iu)(l — Z(ui)) = — z7tT^~~^ X(uj — uj')sgn{w') 

Uj' 

+itujirT 
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The presence of Z{ui) in the denominator of 
C(oj) means that we should redefine A(oj) = 
TXL<lV,a/)<Mw')/ZH in Eq. (§ so that A(w) is 
given by 



Z(w)A(w) = ttT £ ^1 
— /i* — t In 
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(10) 



Equations (^|) and ([10]) constitute the main result of this 
paper. They are the Eliashberg equations for strong cou- 
pling superconductors including the effects of Coulomb 
interaction and disorder. 

The above equations have been derived using the 
pair propagator approach, whereas they are normally de- 
rived in the framework of the Nambu-Gor'kov self-energy 
formalism. We will now extend the latter formalism to 
include the effects of disorder. The various contributions 
to the superconducting self-energy, 



S(w) = iuj[l - Z(u)]t + Z(lo)A(lo)t u 



(11) 




(e) 



FIG. 2. Contributions to the superconducting self-energy 
due to: (a) impurity scattering, (b) electron-phonon interac- 
tion, (c) bare Coulomb repulsion, (d) and (e) disorder-dressed 
Coulomb interaction. 

are shown in Fig. (2). The phonon and bare Coulomb 
contributions are 



^ r ,\ tVw ,Ju't - A(ui')n 



W 
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where W = v /w 2 + A 2 (w) and W = ^/uj' 2 + A 2 {u'). 
They can be found in Ref. |l(|. The disorder-dressed 
Coulomb diagrams of Fig. (2d) and (2e) give 
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In the limit where T — > T Cl the order parameter 
A{u) — » 0, and we obtain the previously derived Eliash- 
berg equations (§) and (|ic|). 

Having derived the Eliashberg equations in Mat- 
subara frequencies, let us now proceed to their numer- 
ical solution. Following Refs. | 10| , ^| , we multiply Eq. 
(§) by A(uj) and equate it to Eq. (jnj. Noting that 



lu = (2n + 1)ttT, J = (2ri + 1)ttT, and defining 
A„ = A(w)/(2n + l) 1 / 2 , we obtain the matrix equation 



where 



^ Snn'A n i — 0, 
n'=0 



(14) 



S nn > = |A(n - n') + X(n + n' + 1) - 2//* - 2t In 
-<y^|l+fA(0)+2^A(n'')Jp r 



M 



n + n' + 1 



[(2n+l)(2n' + l)]V2 



_J ^ 1 ) 

l)(2n' + 1)]V2 + i „ + n » + i 



(15) 
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weak coupling theories is unclear. Our approach has the 
advantage of yielding a compact set of equations that 
are easy to use, and whose weak coupling and weak dis- 
order limits reproduce the known results. In future work 
we will analytically continue these equations to real fre- 
quencies, and use them to analyse the effect of disorder 
on tunneling characteristics in strong coupling materials. 
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FIG. 3. Relative suppression of the transition tempera- 
ture T c vs t ln 3 (l/27rT Cj or) for different values of the elec- 
tron-phonon coupling A. All calculations are performed using 
an Einstein spectrum. 

and M = 1/2ttT c t. The numerical strategy is to find the 
temperature at which the largest eigenvalue of the ma- 
trix S becomes positive. We have to specify the phonon 
elements, A(n), the Coulomb pseudopotential, /i*, the 
effective disorder strength, t = Rb/Rq, and the elastic 
scattering rate, 1/r. As an illustration, we solve the 
equation for an Einstein phonon spectrum, a 2 F(z) — 
\oj e 6(z-uje)/2, for which X(n) = \uj 2 e /[u 2 e + (2mrT) 2 ]. 
We plot T c (A,i)/T c (A,0) vs tin 3 (1/2ttT c (A, 0)t), the lat- 
ter being the initial slope predicted from perturbation 
theory. This enables us to see if the shape of the curve 
deviates from the weak coupling prediction as A is in- 
creased. The only parameters we fix are /i* = 0.1 and 
1/ujet = 250 which arc typical values for lead films. The 
results are shown in Fig. 3. We see that there is a signifi- 
cant deviation from the weak coupling predictions which 
may be important in fitting experimental data. We in- 
tend to extend our numerical analysis to lower values of 
A and different phonon spectra. 

In conclusion we have developed a consistent treat- 
ment of the effect of disorder and strong coupling on the 
transition temperature of superconducting films. Previ- 
ous authors have considered similar systems within an 
exact eigenstate framework [ [THI , but this approach re- 
mains somewhat opaque physically, and its relation to 
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